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ABSTRACT 

We present evidence for the existence of substructure in the "relaxed appearing" cold front clusters 
Abell 2142 and RXJ1720. 1+2638. The detection of these substructures was made possible by com- 
prehensive multi-object optical spectroscopy obtained with the Hectospec and DEep Imaging Multi- 
Object Spectrograph instruments on the 6.5m MMT and 10m Keck II telescope, respectively. These 
observations produced 956 and 400 spectroscopically confirmed cluster members within a projected 
radius of 3Mpc from the centers of Abell 2142 and RXJ1720. 1+2638, respectively. The substructure 
manifests itself as local peaks in the spatial distribution of member galaxies and also as regions of 
localized velocity substructure. For both Abell 2142 and RXJ1720. 1+2638, we identify group-scale 
substructures which, when considering the morphology of the cold fronts and the time since pericen- 
tric passage of a perturber estimated from the cold front radii, could plausibly have perturbed the 
cluster cores and generated the cold fronts observed in Chandra images. The results presented here 
are consistent with cold fronts being the result of merger activity and with cold fronts in relaxed 
appearing clusters being due to minor merger activity. 

Subject headings: galaxies: clusters: individual (Abell 2142, RXJ1720. 1+2638) — X-rays: galaxies: 
clusters 



1. INTRODUCTION 

One of the triumphs of modern cosmology is the suc- 
cess the hierarchical formation model has had in explain- 
ing the observed large scale structure of galaxies. In this 
model, small-scale perturbations in the matter density 
distribution at early epochs intensify due to the effects of 
gravity, collapse and merge together to form the galaxies, 
groups, and clusters observed at later epochs. Since clus- 
ters are the structures to form last, the process of hierar- 
chical formation is expected to reveal itself dramatically 
in such systems. This is borne out by observations, where 
many clusters in the local universe are found to have re- 
cently undergone significant growth via either a major 
merger with another cluster or through minor mergers 
with smaller clusters and groups of galaxies. This offers 
a unique opportunity to study the effects of large scale 
structure formation on the constituent galaxies and intra 
cluster medium (ICM). However, both the detection and 
distinction of clusters undergoing major or minor merg- 
ers can be an observationally expensive exercise in the 
absence of easily de tectable signatures of each process 
(jOwers et al.ll2009bD . 

In this vein, the "cold fronts" first de tected and charac- 
terized by the Chandra X-r ay satellite (|Markevitch et all 
120001 iVikhlinin et al.lfeOOll ). thanks to its excellent sen- 
sitivity and unrivalled spatial resolution, may play an 
important role. Cold fronts are contact discontinu- 
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ities in the ICM which form at abrupt interfaces be- 
tween low and high en t ropy plasma (for a review see 
Markevitch fe V ikhlininl l2007f). In the scenario put 
forth by iMarkevitch et al.l ( 2000T) . commonly referred 
to as the "remnant core" scenario, the edges delin- 
eate the low entropy gas of two cool cores, which 
have just survived a pericentric passage, and the sur- 
rounding hot shocked atmospheres which were stripped 
from the cores during the merger. This remnant core 
scenario is observed in a number of clusters, e.g., 
1ES0657-558 (otherwis e known as the Bullet cluster; 
Markevitch et al.l l2002D. Abell 1758 (IDavid fc Kempnerl 



2004H . Ab ell 115 (IGutierrez fc KrawczvnsklT feoOSL 
Abell 2744 (iKempner fc pavidll2004 lOwers et al.lboilf) . 
Abell 2146 feussell et al.l l2010li and probably Abell 2069 
(|Owers et alj l2009bf ). As the number of clusters ob- 
served by Chandra increased, two things emerged (i)- 
The remnant core scenario cannot explain the majority 
of the observed cold fronts, and (ii)-There is a clear di- 
chotomy when considerin g the X-ray morphology of clus- 
ters hosting cold fronts (IMarkevitch fc V ikhlinin 120071 
lOwers et al1l2009bllGhizzardi et al.ll2010f ). One subset of 
cold front clusters exhibit clearly disturbed X-ray mor- 
phologies, while the second subset have cool cores and 
appear relaxed at large radii. The disturbed morphology 
present in the former subset is clearly related to ma- 
jor merger activity which manifests itself at other wave- 
lengths. The cold fronts due to the remnant core scenario 
are generally found in these clusters, but other merger- 
related processes also produce cold fronts (e.g., the "ram 
pressure slingshot" feature in Abell 168 reported by 
lHallman fc Markevitchl 12004 see also iPoole eFall 120061 
for different cold fronts produced in merger simulations). 
In these cases, cold fronts are excellent signatures of ma- 
jor merger activity (jOwers et al.ll2009a| cT). For the subset 
with relaxed large scale X-ray morphologies, the best ev- 
idence for merger activity comes from the observed cold 
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fronts — the evidence for merger activity at other wave- 
lengths in these clusters is generally ambiguous or non- 
existent. 

The leading hypothesis for the formation of cold fronts 
in clusters with a regular X-ray morphology is the 
relative motion of the low entropy co re gas — termed 
"sloshing'Q by iMarkevitch et aLl ([2001D — with respect 
to the quasi-static clu s ter po t ential. The simulations 
of iTittlev fe HenriksenI (12005T). lAscasibar fc Markevitchl 
(2006ft and lRoediger et aLl ([20111) showed that a gravita- 
tional perturbation in the form of an infalling subcluster 
without a gaseous atmosphere can induce long-lasting os- 
cillations of the the dark matter and gas cores. While the 
precise details differ (e.g., in the lTittlev fc HenriksenI sce- 
nario, t he dark matter oscillates, draggin g the gas with it, 
while in IAscasibar"fe Mark evitchl ([20061 ) the dark matter 
is nearly static while the gas oscillates) in all cases the 
relative motion of low entropy cool core gas with respect 
to the cluster on larger scales produces cold fronts where 
it meets the higher entropy gas at larger radii. Hydro- 
dynamical perturbat ions in the form of m inor merger- 
related weak shocks dChurazov et aL|l2003[) . shocks from 
fast moving galaxies (jRoediger et al.ll2011l ). or acoustic- 
gravity w aves related to bulk motions outside of the clus- 
ter cores (Fuii ta et al.ll2004l ) have also been investigated 
by simulations. In these scenarios, the core gas is dis- 
placed from the static dark matter core and oscillates 
around the potential minimum as the system relaxes back 
to its equilibrium configuration. Several authors have 
suggested that outbursts from AGN m ay provide a suffi- 
cient perturbation to induce sloshing ([Markevitch et al.1 
120011 iHlavacek-Larrondo et al.ll201l( ) although this pos- 
sibility is yet to be fully explored. The emerging con- 
sensus is that the cold fronts seen in relaxed appearing 
cluster are due to minor merger activity. This leads to 
the interesting implication that cold fronts can be used as 
both a signpost of merger activity and, in combination 
with X-ray morphology, a tool to differentiate clusters 
undergoing major and minor mergers. 

Much of the focus on sloshing cold front clusters has 
been at X-ray wavelengths or has involved simulations 
designed to explore mechanisms which are capable of 
reproducing the cold fronts along with a relaxed mor- 
phology at larger radii. To date, no studies have fo- 
cused on identifying possible perturbers by utilizing the 
spatial and kinematic information provided by multiob- 
ject spectroscopy (MOS) — an issue critical to the minor 
merger scenario where a recently merged perturber is 
not apparent at X-ray wavelengths, but may manifest 
itself as a dynamical substructure. Here, we set out 
to address this problem using the powerful combination 
of spatial and redshift information provided by multi- 
object spectroscopy in an attempt to detect merger re- 
lated substructure and characterize its dynamical prop- 
erties. Such observations are a crucial test for constrain- 
ing simulations and for showing that cold fronts are ex- 
cellent probes of ongoing merger activity. In this pa- 
per we present and analyze comprehensive MOS obser- 
vations of the sloshing cold front clusters Abell 2142 
and RXJ1720.1+2638 (hereafter A2142 and RXJ1720, 

4 For clarity, throughout the paper we will refer to clusters which 
have regular/relaxed X-ray morphologies with cold fronts in their 
cores as sloshing cold front clusters. 



respectively) which were selected from th e subsample 
of relax ed appearing cold front clusters in lOwers et al.1 
(|2009bD . Previous MOS follow up of cold front clusters 
taken from the disturbed subset in lOwers et al.l ([2009b) 
has demonstrated that MOS observations are an essen- 
tial ingredient requ ired for understand ing the merging 
history of a cluster (jOwers et al.ll2009"cT ia!). 

A2142 was one of the first clusters targeted with the 
Chandra X-ray telesco pe and, along with t he obser- 
vations of Abell 3667 (|Vikhlinin et al.1 120011 ). charac- 
terized for the first tim e the cold front phenomenon 
([Markevitch et ail I2000D. The Einst ei n X-r ay images 
of A2142 presented in lOegerle et al.1 JT995), revealed 
an asymmetric X-ray morphology, while their analy- 
sis of 103 spectroscopically confirmed members of this 
cluster yielded a mean redshift of z = 0.0905 ± 
0.0004 and a velocit y dispersion a = 1280tygkms _1 . 
lOegerle et al.1 (119951) performed a Dressler-Schectman 
test ([Dressier fc Shectmanl I1988D and found marginal 
evidence for dynamical substructure, while their data 
also revealed a high peculiar velocity for the second 
ranked member of ~ 1600 km s . X-ray observations 
using ROSAT revealed a highly elliptical X-ray mor- 
phology with a luminous compact core, the centroid of 
which is offset from that of the larger s cale emission 
([Buote fc Tsal Il996l IHenrv fc Briell pL996ft . The tem- 
perature map of IHenrv fc Briell ([1996 ) showed that this 
compact core was cooler than the surrounding gas and 
that the temperature distribution was asymmetric on 
large scales. The high ellipticity, offset core, asymmet- 
ric temperature distribution and evi dence for dynam ical 
structure wer e interpreted by both IHenrv fc Briell and 
iBuote fc Tsail as evidence for a later stage merger (i.e., 
viewed ~ 1 Gyr after the merger) in A2142. Hints of 
a substructure to the northwest of the core in A2142 
ar e seen in the weak l ensing and SZ maps presented 
in lUmetsu et al.[ (12009). and also the lensing data of 
Oka be fc Umetsul ([20081 ). 

RXJ1720 was categorized as an extended X-ray source 
in the Einstein Slew Survey ([Elvis et al.1 119921 ) . The 
Chandra observations of RXJ1720 revealed two cold 
fronts within 100" of the cluster center, while the X-ray 
emission outside of this radiu s showed a regular mor- 
phology (Maz zotta et al.lfeOOll) . These Chandra observa- 
tions of RXJ1720 provided the first clear indication that 
cold fronts c an also be found in re laxed appearing cool 
core clusters. iMazzotta et al.l ([20011 ) interpreted the cold 
fronts as evidence for either a late stage major merger 
viewed after several pericentric passages or the asyn- 
chronous collapse of two nearly cospatial density pertur- 
bations with diffe rent linear scales. Using a deepe r Chan- 
dra observation, IMazzotta fc Giacintuccil (|2008| ) found 
that the spiral features revealed in the temperature and 
residual maps are more readily explai ned by the slosh- 
ing scenario. [Mazzotta & Giacintuccil also reported an 
interesting correlation between the spiral structure and 
radio emission from a mini-halo and suggested that the 
origin of this diffuse radio emission may be related to 
the sloshing process. At optical wavelengths, RXJ1720 
was detected with the C4 cluster finding a lgorithm us- 
ing t he Sloan Digital Sky Survey (SDSS) ([Miller et al.1 
120051) . From 26 sp e ctrosc opically confirmed cluster mem- 
bers I Miller et al.1 ([20051 ) measured a mean redshift of 
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z = 0.1603 and a velocity dispersion a — 878 kms" 1 
while the Dressier- Schectman test revealed no evidence 
for dynamical substructure. The p rojected mas s maps 
from the weak lensing analysis of Ok abe et al.l ()2010ft 
showed that RXJ1720 harbors a second mass conden- 
sation just north of the main cluster component which is 
indicative of merger activity. 

For both clusters there are tantalizing hints for merger 
activity. The goal of this paper is to seek conclusive evi- 
dence of dynamical substructure related to merger activ- 
ity by using the extensive spectroscopy collected at the 
6.5m MMT and 10m Keck II Telescope using the Hec- 
tospec and DEep Imaging Multi-Object Spectrograph in- 
struments, respectively. We describe our target selection 
and MOS observations in Section [2] and our methods for 
cluster member selection and substructure detection in 
Section [3] We interpret our results in Section 2] and 
present a summary and concluding remarks in Section [5j 
Throughout this paper we assume a standard ACDM 
cosmology where Hq = 70 kms -1 Mpc -1 , f2 m = 0.3 and 
Oy = 0.7. Assuming this cosmology, at the systemic red- 
shifts of A2142 and RXJ1720, 1 arcsecond corresponds 
to a physical scale of 1.68 kpc and 2.76 kpc, respectively. 

2. OBSERVATIONS AND DATA REDUCTION 

2.1. Photometric Catalogs 

The spectroscopic target catalogs for both A2142 and 
RXJ1720 were drawn from the SDSS and included ob- 
jects within a projected distance from the cluster cen- 
ter of i? proj = 3Mpc, or - 30' and - 18' for A2142 
and RXJ1720, respectively. The catalogs were fil- 
tered to contain only objects identified as galaxies by 
the SDSS pipeline mo rphological classification scheme 
flStoughton et all 12002ft . Objects near to bright stars 
and misclassified as galaxies by the automated SDSS 
pipeline were identified by eye and rejected from the 
catalog. Objects with existing SDSS spectroscopy and 
having redshifts which placed them well outside the red- 
shift range expected for the cluster of interest were also 
removed from the sample. We transformed the SDSS 
g and r— band magnitudes onto the Johnso n-Cousins B 
and R system using Equations A5 and A7 of lCross et all 
(l200l . 

2.2. MMT/Hectospec Observations 

The majority of the spectra used in this paper were 
collected at the 6.5m MMT using the Hect ospec multi- 
object spectrograph (|Fabricant et all [2005ft . Hectospec 
is a bench-mounted spectrograph which has the capabil- 
ity to simultaneously observe 300 objects through 1".5 
diameter fibers spread over a 1 degree field of view. The 
observations were performed with the 270 groove mm -1 
grating resulting in spectra with ~ 6 A resolution cover- 
ing the wavelength range 3500— 10000 A. Approximately 
30 fibers were allocated to blank sky regions to facilitate 
sky subtraction during the data reduction. Due to the 
large number of targets compared to the number of avail- 
able fibers and fiber packing limitations, 6 and 8 config- 
urations were required to obtain adequate spectroscopic 
completeness for RXJ1720 and A2142, respectively. The 
exposure time, magnitude limit and seeing conditions for 
each MMT/Hectospec configuration are summarized in 
Table [1] 



The fiber configurations were generated with the 
XFITFIBS0 software where the ranking capability was 
utilized to ensure that the most likely cluster members 
had a better chance of being allocated to a configura- 
tion. For RXJ1720, objects were ranked based on their 
position on the color-magnitude relation (Figure [lj and 
also their i? pro j in the sense that objects lying redward 
of the cluster red-sequence, i.e., above the bold black 
line in Figure [I] and with 2.5 < -R pro j < 3Mpc were 
ranked lowest while those galaxies which lie on or blue- 
ward of the red-sequence and with -R pro j < 0.5 Mpc were 
ranked highest. The slope of the line delineating those 
galaxies redward of the red sequence is calculated from 
the best fit to the red s equence slope versus reds hift di- 
agram in Figure 3. of lLopez-Cruz et all ([2004ft . The 
ranking procedure was similar for the A2142 configura- 
tions, except that objects lying redward of the cluster 
red-sequence were simply not included for observation 
(Figure [J), while a number of objects identified as be- 
ing part of local projected overdensities (aside from the 
central overdensity) were given higher rankings and a 
number of fainter (20 < R < 20.5) objects were included 
with low rankings as filler objects to be allocated to spare 
fibers. 

The spectra were reduced at the Smithsonian Astro- 
physical Observatory Telescope Data_Centei0 (TDC) us- 
ing the SPECROAD pipelin£] (jMink et all (2007ft . Red- 
shifts were also determined at the T DC using the IRAF 
cross-correlation XCSAO software (K urtz et all 11992ft 
and each spectrum was assigned a redshift quality of "Q" 
for a reliable redshift, "?" for questionable and "X" for a 
bad redshift measurement. These observations provided 
1771 and 815 quality "Q" redshift measurements for ex- 
tragalactic objects in the A2142 and RXJ1720 fields, re- 
spectively. 

2.3. Keck/DEIMOS Observations 

The MMT/Hectospec data for RXJ1720 were sup- 
plemented with spectra obtained at the Keck II tele- 
scope using the Deep Imaging Multi-Object Spectro- 
graph (DEIMOS). In MOS mode, DEIMOS utilizes a 
slitmask to collect spectra for ~ 100 target objects po- 
sitioned within a ~ 16.7' x 5' field of view. The observa- 
tions presented here made use of the 600 line mm -1 grat- 
ing which, in combination with a 1".5 slitwidth, gives a 
spectral resolution of ~ 7.0 A. The GG400 order-blocking 
filter was used and the grating was tilted to a central 
wavelength of 5500 A which yields a useable spectral 
range of - 4000 - 9000 A. 

Slit masks were generated with the DSIMULATOR 
task within the DEIMOS IRAF package. The input tar- 
get catalog was limited to contain objects with R < 22 
and within i? pro j= 2 Mpc (~ 12'). The catalog was culled 
of those objects found to be non-cluster members from 
the previous MMT/Hectospec observations, while those 
spectroscopically confirmed cluster members with high 
quality Hectospec spectra (i.e. S/N> 8) were also re- 
moved. The DSIMULATOR software allows a two tier 
ranking of targets: First, the catalog is split into subsam- 
ples with objects in the first subsample given first priority 

5 http:/ /www. cfa.harvard.edu/mmti/hectospoc. html 

6 http://tdc-www.harvard.edu 

7 http:/ /tdc-www. harvard.edu/instruments/hectospec/specroad. html 
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Figure 1. Color magnitude diagrams for the objects within 3Mpc of A2142 (left panel) and RXJ1720 (right panel). Filled green circles 
are spectroscopically confirmed cluster members, blue stars are foreground objects, red crosses are background objects and black crosses 
are objects without a spectrum. The black lines show the estimates used for t he re d and blue limits defining the cluster red envelope which 
was used f or th e cl uster membership allocation procedure described in Section l3.f I The red limit was used in the ranking process described 
in Sections [2?2l and [231 



Table 1 

Summary of the MMT/Hectospec observations. 



Cluster 


Date 


Magnitude limit 


Exposure Time (s) 


Nspcc 


N z 


Seeing 


RXJf720 


2006 May 4 


R < 20.5 


3 X 900s 


208 


189 






2006 May 4 


R < 20.5 


3 x 900s 


186 


176 






2007 Mar 15 


R < 20.5 


3 x 900s 


189 


171 


0".97 




2007 Mar 16 


R < 20.5 


3 x 1200s 


177 


92 


0".63 




2007 April 16 


R < 20.5 


3 x 1200s 


161 


122 


1".5 




2007 April 20 


R < 20.5 


5 x 600s 


166 


121 


1".3 


A2142 


2009 February 5 


R < 20.5 


3 x 1200s 


266 


255 


1".08 




2009 February 21 


R < 20.5 


3 x 1200s 


261 


174 


1".30 




2009 February 24 


R < 20.5 


1 x 1200s 


263 


159 


0".56 




2009 March 1 


R < 20.5 


3 x 1200s 


263 


249 


0".84 




2009 April 22 


R < 20.5 


3 x 960s 


260 


255 


0".60 




2009 April 23 


R < 20.5 


3 x 1200s 


259 


244 


1".21 




2009 April 23 


R < 20.5 


3 x 1200s 


255 


240 


0".61 




2009 April 24 


R < 20.5 


3 x 1200s 


212 


229 


0".86 



in slit allocation. Second, the objects within subsamples 
are given weights such that the highest weighted objects 
are given a higher priority during slit allocation. We di- 
vide the catalog into 5 subsamples: the hrst and second 
contain objects lying on or blueward of the red sequence 
with R < 21.5 and 21.5 < R < 22, respectively, the third 
contains objects already observed with Hectospec, the 
fourth and fifth contain objects lying redward of the red 
sequence with R < 21.5 and 21.5 < R < 22, respectively. 
The object weights were assigned according to i? pro j- 
This subsampling and weighting scheme means that ob- 
jects within the central 500 kpc lying on or blueward of 
the red sequence and with R < 21.5 are most likely to 
be allocated a slit, while objects with 21.5 < R < 22 
lying redward of the red sequence and 2 Mpc from the 
cluster center are least likely to be allocated a slit on the 
mask. The position angle of each slit may be tilted by 
±30 deg with respect to the mask position angle (which 
is aligned with the spatial direction on the CCD). Where 
the position angle of a galaxy's major axis (taken from 
the SDSS shape parameters) does not exceed this limi- 



tation, we tilt the slit to align with the major axis. In 
practice, this was only possible in a few cases and the 
majority of slits were therefore assigned a position angle 
of 5 deg with respect to the mask. Using this setup, we 
generated 6 slitmasks which were observed on 2009 May 
24. Table [5] summarizes these observations. 

The data were reduced with the spec2d software^ devel- 
oped for the DEEP2 galaxy redshift survey ()Davis et al.l 
2003). BriefljQ the software uses the internal Quartz 
flats taken on the afternoon before the observations to 
find and rectify the slits on the CCD, to determine the 
throughput of each slit and to correct for fringing. We 
utilize blue and red specific arc lamps to determine the 
2D wavelength solution which is initially estimated using 
the DEIMOS optical model. The science exposures are 
combined using an inverse variance weighed mean after 
cosmic ray correction and sky subtraction. The I D spec- 
tra are then extracted using optimal extraction (jHornd 

8 http:/ /deep. berkeley.edu/spec2d/ 

9 a more thorough explanation ca n be found at 
http: / /deep. berkeley.edu/spec2d/primer. html 
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Table 2 

Summary of the Keck/DEIMOS observations of RXJ1720. 



a (J2000) 


<5 (J2000) 


Mask PA 


Exposure Time (s) 


Nspcc 


N z 


Seeing 


17:20:05.02 


26:39:41.5 


180.0 


2 X 1100 -4 


- 1340s 


70 


61 


0" 


.57 


17:20:14.28 


26:35:38.2 


0.0 


2 x 1100 4 


- 1340s 


67 


67 


0" 


.60 


17:20:03.70 


26:37:31.9 


0.0 


2 x 1100 4 


- 1200s 


62 


57 


0" 


.50 


17:20:16.73 


26:37:25.3 


180.0 


2 x 1100 4 


- 1200s 


60 


57 


0" 


.63 


17:19:46.02 


26:37:31.9 


0.0 


2 x 1100 4 


- 1500s 


68 


68 


0" 


.54 


17:20:34.44 


26:37:31.9 


180.0 


2 x 1100 - 


4- 800s 


62 


50 


0" 


.61 



1986) and rebinned to linear wavelength bins using the 
IRAF task RSPECTEXT. 

The redshifts are again measured with the XCSAO 
cross-correlation package in IRAF (K urtz et al.l 119921) 
and each spectrum was visually assessed (by MSO) and 
assigned a redshift quality of "Q" for a reliable redshift, 
"?" for questionable and "X" for a bad redshift mea- 
surement. These observations produced 333 quality "Q" 
redshift measurements for extragalactic objects in the 
RXJ1720 field. 

2.4. Redshift precision 

The redshift precision is determined here using repeat 
observations taken during the observating runs. For the 
MMT/Hectospec observations of A2142, there were 234 
objects having two quality "Q" redshift determinations 
which were reobserved due to their lower signal to noise 
spectra (none were taken for RXJ1720). Taking the dif- 
ference of the redshifts of these repeat observations (in 
the sense that we subtract the subsequent higher S/N 
redshift measurement from the lower S/N initial red- 
shift) we use biweight estimators to find a mean redshift 
difference of < Acz >— 2.3 ± 2.3 kms -1 and a stan- 
dard deviation a(Acz) = 34.2 ± 2.3 kms -1 indicating a 
single measurement uncertainty of c(< Acz >)/\/2 = 
24.2 ± 1.6 kms -1 . This single measurement uncertainty 
is consistent with the median of the redshift uncertainties 
determined by the XCSAO software of 24.7 kms -1 . 

We can test the veracity of this internal precision mea- 
surement by comparing the repeat observations taken 
with MMT/Hectospec and from the SDSS. There are 
301 objects which have redshift measurements from both 
our A2142 MMT/Hectospec observations and the SDSS. 
Taking the difference between the MMT/Hectospec and 
SDSS redshifts, we find < Acz >= -4.8 ± 1.8 km s -1 
and a(Acz) = 31.3 ± 1.5 kms -1 . This standard devi- 
ation is well encompassed by the median of the indi- 
vidual redshift uncertainty measurements for the SDSS 
(czcrr = 47.5 kms -1 ) and the MMT/Hectospec (cz er r = 
21.9 kms -1 ) observations. Similarly, while there were 
no repeat observations within the MMT/Hectospec or 
Keck/DEIMOS redshift catalogues for RXJ1720, there 
were 52 common measurements between the SDSS and 
MMT/Hectospec. The mean difference in the redshift 
measurements common to the SDSS and Hectospec cat- 
alogues was < Acz >= — 31. 6± 7.4 kms -1 and a(Acz) = 
49.8 ± 6.3kms -1 . The quadrature sum of the me- 
dian value for the individual uncertainties on the SDSS 
(cz orr = 45.8 kms -1 ) and the MMT/Hectospec (cz eT i = 
28.9 kms -1 ) repeat observations well encompasses the 
scatter in < Acz > implying that the individual uncer- 



tainties provide a good measure of the real uncertainties 
associated with these measurements. 

Given the robustness of the MMT/Hectospec redshift 
measurements and associated uncertainties, we now use 
the 33 MMT/Hectospec objects which were reobserved 
with Keck/DEIMOS to determine the real uncertain- 
ties in the DEIMOS measurements, finding < Acz > — 
-69.5 ± 19.0kms -1 with scatter a(Acz) = 111.6 ± 
18.5 kms -1 . Comparing the median values of the indi- 
vidual uncertainty measurements (for MMT/Hectospec: 
cz^7- = 38.5 kms -1 and for Keck/DEIMOS: cz^ = 
52.3 kms -1 ) it can be seen that the scatter is not well 
accounted for by the quadrature sum of the uncertain- 
ties. We conclude that the real redshift uncertainties on 
the DEIMOS measurements are ~ 100 kms -1 . 

Since there are significant offsets measured in the re- 
peat observations for RXJ1720, when combining the 
MMT/Hectospec, Keck/DEIMOS and SDSS catalogs we 
correct each Keck/DEIMOS and SDSS redshift measure- 
ment by the < Acz > value measured above such that 
they are consistent with the MMT/Hectospec measure- 
ments. Where there are repeat measurements, for both 
the A2142 and RXJ1720 catalogs we use the measure- 
ment with the lowest measured uncertainty. Where a 
Keck/DEIMOS measurement is involved, we double the 
redshift uncertainty, consistent with the results outlined 
above. The final catalogs contain 1121 and 1635 single 
redshift measurements for extragalactic objects in the 
field surrounding RXJ1720 and A2142, respectively. We 
tabulate the position, redshift, redshift uncertainty and 
redshift source of objects in the final catalogs in Tables [3] 
and H for A2142 and RXJ1720, respectively. 

2.5. Spectroscopic completeness 

An understanding the spectroscopic completeness of 
our redshift catalogs is imperative when interpreting the 
results of tests for substructure, in particular for those 
that rely on finding local overdensities in the spatial 
distribution of the cluster members. More specifically, 
we need to understand what fraction of potential clus- 
ter members are in the photometric catalog but have 
not been assigned a redshift. Therefore, we define the 
spectroscopic completeness as the ratio of the number of 
objects with quality "Q" redshift measurements to the 
number of objects in the photometric catalog. As seen 
in the right panel of Figure [TJ objects lying red ward of the 
cluster red sequence are unlikely to be cluster members. 
Therefore, given our goal of understanding how redshift 
incompleteness affects the cluster members, for RXJ1720 
we only measure the spectroscopic completeness for the 
subsample of objects lying below the upper black bold 
line defining the red limit of the cluster red sequence in 
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Table 3 

Combined redshift catalog for A2142. 



R.A. (J2000) decl. (J2000) (cz) cz uncertainty cz source 



deg. 


deg. 


km s 1 


kms 1 




239.036453 


27.235901 


57055.78 


18.59 


MMT/Hectospec 


239.042587 


27.177689 


28939.35 


30.64 


MMT/Hectospec 


239.045349 


27.174765 


28811.72 


56.53 


MMT/Hectospec 


239.048401 


27.247974 


96294.95 


58.26 


MMT/Hectospec 


239.049500 


27.187420 


24531.00 


48.49 


MMT/Hectospec 


239.053802 


27.157728 


86498.82 


35.92 


MMT/Hectospec 


239.055679 


27.257549 


33359.37 


35.87 


MMT/Hectospec 


239.062576 


27.132755 


97387.20 


21.20 


MMT/Hectospec 


239.078674 


27.185207 


67723.70 


27.55 


MMT/Hectospec 


239.079071 


27.307705 


26771.70 


19.48 


MMT/Hectospec 



Note. — This tabic is available in its entirety in a machine-readable form in the 
online journal. A portion is shown here for guidance regarding its form and content. 



Table 4 

Combined redshift catalog for RXJ1720. 



R.A. (J2000) decl. (J2000) (cz) cz uncertainty cz source 



deg. 


deg. 


kms 


km s 1 




259.713308 


26.620569 


119987.97 


102.64 


MMT/Hectospec 


259.718571 


26.644934 


47540.42 


144.19 


MMT/Hectospec 


259.719717 


26.702864 


111920.71 


101.03 


MMT/Hectospec 


259.726667 


26.565548 


72713.53 


28.23 


MMT/Hectospec 


259.727325 


26.544941 


97930.56 


28.17 


MMT/Hectospec 


259.728700 


26.560825 


72656.78 


80.44 


MMT/Hectospec 


259.731187 


26.554123 


98086.70 


55.13 


MMT/Hectospec 


259.732817 


26.646864 


36395.72 


35.47 


MMT/Hectospec 


259.738025 


26.516831 


30848.81 


21.77 


MMT/Hectospec 


259.739112 


26.738709 


117214.30 


35.72 


MMT/Hectospec 



Note. — This table is available in its entirety in a machine-readable form in the 
online journal. A portion is shown here for guidance regarding its form and content. 



the right panel of Figure[T] This also allows a more mean- 
ingful comparison to the A2142 data. Figure [5] shows the 
spectroscopic completeness as a function of i? pro j for both 
clusters for different magnitude bins. For A2142, it can 
be seen in the left panel of Figure [2] that for the high pri- 
ority targets (i.e. those with R < 20.) the spectroscopic 
completeness is always higher than 80%, while it drops to 
30 — 40% for the low priority faint objects. We therefore 
achieve a high level of spectroscopic completeness to ~ 3 
magnitudes fainter than the characteristic magnitude of 
the cluster luminosity function, M*. 

With regard to RXJ1720, the right panel in Figure [2] 
shows that we achieve spectroscopic completeness above 
~ 85% for magnitudes brighter than R = 20.5 (apart 
from magnitudes 19.5 < R < 20.5 at radii outside 
2Mpc). At fainter magnitudes than R = 20.5 our data 
are limited to the central ~ 2 Mpc region covered by the 
Keck/DEIMOS spectroscopy. For these observations, the 
smaller FOV and more limited number of objects ob- 
served per pointing (compared to the Hectospec obser- 
vations) meant that it was difficult to obtain multiple 
pointings of the same regions over the entire central area 
within 2 Mpc of the cluster center. Therefore, the spec- 
troscopic completeness at these faint magnitudes is lim- 
ited, but still respectable, at > 50% for radii < 1.5 Mpc 
and 20.5 < R < 21.5. As expected, the low prior- 
ity 21.5 < R < 22.0 objects have low spectroscopic 
completeness (~ 40% for radii < 1.5 Mpc). Given the 



more heterogeneous coverage for RXJ1720, we further 
explore the spatial distribution of the spectroscopic com- 
ple teness. Briefly, we use the WVT binning algorithm 
bv iDiehl fc Statlerl (|2006[). which is a generalization of 
Cappellar i fc Copinl 's (|2003l ) Voronoi binning algorithm, 
to produce an adaptively binned image of the spatial dis- 
tribution of objects in the photometric subsample such 
that each bin contains at least 20 galaxies. This binning 
is then applied to the spatial distribution of objects in 
the spectroscopic sample. The resulting image is divided 
by the binned image of the spatial distribution of objects 
in the photometric sample to produce the completeness 
map. These maps are produced at three limiting mag- 
nitudes of interest: R < 20.5, R < 21.5 andi? < 22.0 
and are shown in Figure [3] The leftmost map in Fig- 
ure [3] shows the completeness for the magnitude limit 
of the MMT/Hectospec observations and confirms that 
the spectroscopic completeness is excellent and generally 
above 80%. The middle and right panels reveal a more 
patchy spectroscopic completeness which is high in the 
central 1 Mpc and falls off towards the cluster outskirts. 

3. ANALYSIS 

3.1. Allocation of cluster members 

We allocate cluster membership using a two-step pro- 
cedure. In the first step, crude cluster membership is 
achieved by selecting galaxies within 3 Mpc of the bright- 
est cluster galaxy (BCG) and with a peculiar velocity of 
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Figure 2. Spectroscopic completeness as a function of radius in different magnitude bins (listed in the top right corner) for A2142 (left 
panel) and RXJ1720 (right panel). 
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Figure 3. Spectroscopic completeness maps for RXJ1720 for objects with magnitudes R < 20.5 (left panel), R < 21.5 (middle panel) and 
-R < 22 (right panel). The dashed circles have radii of 1 and 2 Mpc and the cross marks the position of the central BCG. 



±10000 km s" 1 with respect to the BCG r edshift. The 
shifting Rapper method used previously in lOwers et al.l 
(|2009al lcf) does not adequately remove the remaining in- 
terlopers. This is because the depth and completeness of 
the A2142 redshift catalog means gaps in velocity phase 
space are more likely to be "filled" , while the proximity of 
a nearby structure in velocity space for RXJ1720 causes 
the shifting gapper method to fail (Figure 0}. There- 
fore, we employ the pa r ametr ic technique introduced by 
Iden Hartog fc Katgertl (|1996f ) where the maximum pos- 
sible velocity at the projected radius of each galaxy is 
calculated given the mass profile of the cluster, which is 
der ived from the virial theorem (for a detailed discussion 
see lWoitak fc Lokas|[2007llWoitak et al.ll2007D . Determi- 
nation of the virial mass relies critically on the assump- 
tion that the particles tracing the underlying gravita- 
tional potential are virialized and that their distribution 
is spherically symmetric. The most likely galaxies for 
which these assumptions apply are elliptical and SO type 
galaxies, while spiral galaxies are more likely to be either 
line of sight interlopers, or an infalling po pulation on ra- 
dial orbits and, therefore, not virialized (B iviano et al.l 
119971 iBiviano fc Katgertl [20041 ). For these reasons, we 
use only galaxies on the red sequence shown in Figure [TJ 



which we assume is primarily composed of elliptical and 
SO galaxies, to determine the mass profiles used in defin- 
ing the velocity limits for cluster membership. Further- 
more, to reduce any biases in the measurements of the 
mean projected separations used for the virial mass es- 
timator, we restrict the red sequence galaxies to a mag- 
nitude range with high spectroscopic completeness (i.e., 
R < 19.5 for RXJ1720 and R < 19.0 for A2142). The 
results of the membership allocation are presented in Fig- 
ure [4] where it can be seen that this method of interloper 
removal does an excellent job. For reference, we high- 
light the galaxies lying on the blue cloud using blue stars 
in Figure |H We define blue cloud galaxies as those ly- 
ing blueward of the lower limit defining the red sequence 
in Figure [TJ Inspection of the distribution of blue cloud 
galaxies confirms that they have a very different phase- 
space distribution when compared to the red sequence 
galaxies, thus justifying our decision not to use them 
when estimating the virial mass. We note that our pro- 
cedure rejects a number of galaxies which lie close to 
the boundaries defined by the maximum velocity pro- 
file but also appear to be well separated from the more 
clear-cut non-members. These consist primarily of blue 
cloud galaxies and, given their large peculiar velocities, 
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are likely to be galaxies on radial orbits on their first pas- 
sage onto the cluster. We prefer the more conservative 
cut here, since these galaxies appear to be distributed 
evenly in phase space and are thus unlikely to be part 
of the substructures we seek to detect, while their inclu- 
sion in the sample may affect the tests we plan to use for 
detecting substructure in the later sections. 

The final cluster member catalogs contain 956 and 400 
members for A2142 and RXJ1720, respectively. From the 
ensemble of memb er redshifts we use biweight estimators 
(jBeers et al.lll990T) to measure cluster redshifts of z c i us = 
0.09005 ± 0.00012 and z dus = 0.16010 ± 0.00018 and pe- 
culiar velocity dispersions oVpcc = 995 ± 21kms _1 and 
cr vpcc = 882 ±29 km s" 1 for A2142 and RXJ1720, respec- 
tively, where the uncertainties are la and a re estimated 
using the jack-knife resampling approach (|Beers et al.l 
[199(1 . 

3.2. Substructure Detection 

Ultimately, our goal in this study is to test the hypoth- 
esis that cold fronts in clusters with otherwise regular 
X-ray morphologies can be caused by the gravitational 
perturb ation of an infalling substructure during a minor 
merger (IMa rkevitc h et al.ll200lLlA scasibar & Mar kevitchl 
20061 IMarkevitdi fc Vikhlininl 120071 IQwers et al.[ ]2009b. 
Roediger et al.l 1201 1H ■ To that end, we use our samples 
of spectroscopically confirmed cluster members to search 
for merger-related substructures in the ID velocity, 2D 
spatial and 3D velocity-plus-spatial distributions. 

3.2.1. Testing the shape of the peculiar velocity distribution 

A first order test of the dynamical regularity of a clus- 
ter comes from the shape of the peculiar velocity distribu- 
tion which, in a dynamically relaxed cluster, is expected 
to take an approximately Gaussian shape. Gross depar- 
tures from a Gaussian shape, particularly in the form of 
a skewness or bimod ality, often indicate th e existence of 
merg er activity (e.g.. iZabludoff et al.lll993l IQwers et al.l 
120111 ). We ch oose to use the Gauss -Hermit e reconstruc- 
tion m ethod (fZabludoff et al.lll993l : see also IQwers et al.l 
|2009cD to test for departures from a Gaussian shape. 
Briefly, the velocity distribution is described by a series 
of Gauss-Hermite functions with the Gauss-Hermite mo- 
ments ho ~ 1 multiplying the zeroth order term, which 
is the best fitting Gaussian with mean V and standard 
deviation S, while the h 3 and /14 terms approximate 
asymmetric and symmetric deviations from a Gaussian. 
Radial changes in the velocity distribution profile, e.g., 
the decrease in the veloci ty dispersion with rad i us ob - 
served in many clusters ()den Hartog fe Katgert 19961) . 
can cause significant /14 terms which are not related to 
merger activity. Our aim is to detect merger related sub- 
structure rather than non-Gaussian shapes caused by a 
changing velocity dispersion profile so, in addition to 
measuring the Gauss-Hermite distribution for the en- 
tire cluster member sample, we also repeat the Gauss- 
Hermite reconstruction for members with 1500 <i? pro j< 
3000,750 <i?p ro j< 1500 and i? proj < 750 kpc. To de- 
termine the significance of the terms, we gener- 
ate 10,000 Gaussian random distributions with the same 
number of data points, mean and standard deviation as 
the observed velocity distribution and measure the 
value for each random velocity distribution. We define 



-P[^3,4] = 2 min{P[/i 3i 4(sim) < /i 3j 4(obs), P[h 3 ^(sim) > 
/i3,4(obs)]} where P[/i3,4(sim) < h 3 ^(obs)] and 
P[/i3 i 4(sim) > /i3.4(obs)] are determined from the num- 
ber of times a ft.34 value smaller and larger than the 
observed term is measured in the 10,000 random 
distributions. The results of this analysis are shown in 
Figure [5] for A2142 and Figure © for RXJ1720. 

The A2142 velocity distributions for the radial ranges 
i? P roj< 3000,1500 <i? proj < 3000 kpc have significant 
positive /13 terms. This indicates that the peak of 
the distribution is more negative and that the distri- 
bution has more power in the positive peculiar veloc- 
ity tails than expected from the best-fitting Gaussian. 
The 1500 <-R pro j< 3000 kpc velocity distribution also 
has a significant negative /14 term, which indicates the 
velocity distribution has a broader, more boxy shape 
than the best fitting Gaussian. The velocity distribu- 
tion in the range 750 <-R pr oj< 1500 kpc does not ap- 
pear to have a shape which differs significantly from a 
Gaussian one, but we note that the measured velocity 
dispersion, S — 1184 km s -1 is larger than that mea- 
sured at other radii. The peculiar velocity distribution 
in the central 750 kpc returns marginally significant neg- 
ative /13 and positive /14 values indicating that there are 
symmetric and asymmetric distortions present, although 
the marginal significances do not allow strong conclusions 
to be drawn. 

The velocity distribution for the entire RXJ1720 mem- 
ber sample (top left panel in Figure [6]) is well approxi- 
mated by a Gaussian — small deviations are present at the 
5 — 6% level in both the ft. 3 and /14 terms but they are only 
significant at the 10% level. The velocity distribution for 
the cluster members with 1500 <i? pro j< 3000 kpc, (top 
right panel in Figure [6|) has a negative h 3 term which 
is significant at the 2% confidence level. This indicates 
that the velocity distribution in this radial bin has a peak 
which is shifted towards positive peculiar velocities and 
a tail with slightly more power at negative peculiar ve- 
locities. The remaining two velocity distributions (lower 
panels in Figure [6]) do not harbor any significant devia- 
tions from a Gaussian shape. 

3.2.2. Substructure in the spatial distribution of galaxies 

Previous studies have shown that clusters in the pro- 
cess of merging may contain substructures which have 
little to no significant detectable impact on the shape of 
the peculiar velocity distributions, particularly in cases 
where t he merger plane is perpendicular to ou r line 
of sight (jPinknev et all 119961 IQwers et al]|2009dR ). In 
these cases, the substructure may be conspicuous as lo- 
cal peaks in the projected surface density of the galaxies 
which indicate the existence of compact, bound subclus- 
ters. Here, we utilize our comprehensive spectroscop- 
ically confirmed cluster member samples to search for 
spatially compact substructures. 

Clusters exhibit a large dynamic range in their pro- 
jected galaxy surface density distributions. Therefore, in 
order to both resolve spatially close substructures and 
minimize the effects of noise in low density regions, an 
adaptive smoothing kernel is required. We use a Gaus- 
sian kernel with FWHM equal to the radius of the Nth. 
nearest neighbor in projection, where N — W (n mem ) 
and n mem is the number of members in the sample, to 
smooth the projected galaxy distribution. The results 
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Figure 4. Phase space diagrams used for defining cluster membership for A2142 (left panel) and RXJ1720 (right panel). The red curves 
show the cluster limits in phase-space and are determined from the cluster mass profile (see the text for a description). Black crosses 
show cluster members which lie on the red sequence, while blue stars within the phase space limits show those cluster members which lie 
blueward of the red sequence. Red boxes show fore- and background galaxies. Red boxes filled with blue stars show non-members which 
are bluer than the cluster red sequence. 
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Figure 5. The velocity distribution for all 956 cluster members (top left panel), 463 members with 1500 <R pro j< 3000 kpc (top right), 
270 members with 750 <i? pro j< 1500 kpc (bottom left) and 223 members with R pro j< 750 kpc (bottom right). In each panel, the best 
fitting Gaussian parameters, the Gauss-Hermite terms h^ and h^, representing the asymmetric and symmetric deviations from a Gaussian 
shape, and their associated level of significance (see text) are given in the upper right corner. The histograms show the observed velocity 
distribution and the bin size is set to 3S/(A r ga i)' 1 '' 3 ^ , where -/V ga i is the number of galaxies listed above. The solid black line shows the 
best fitting Gaussian, while the dashed red line shows the Gauss-Hermite reconstriction of the velocity distribution. 
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Figure 6. The same as Figure[5]but for RXJ1720. The number of galaxies in each subsample shown in the top left, top right, bottom left 
and bottom right are 400, 146, 117 and 137 galaxies, respectively. 
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of the smoothing are shown in the top left panels of 
Figures and M for A2142 and RXJ1720, respectively. 
Both clusters exhibit a dense core of galaxies surround- 
ing the central BCG, elongated morphologies on larger 
scales and a number of local overdensities. The central 
regions of A2142 appear to be irregular, with a number of 
local overdensities within 1 Mpc of the cluster center, the 
most significant of which occurs ~ 180 kpc to the NW, 
roughly coincident with the second BCG. There are also 
four large conglomerations of galaxies in the outskirts of 
A2142; one ~ 2.4 Mpc to the SE, one ~ 1.7 Mpc due 
south, one ~ 2.5 Mpc to the SW and one 2.2 Mpc to the 
NW. Most notably for RXJ1720, there is a second peak 
in the projected galaxy density slightly east of due north 
and ~ 550 kpc from the cluster center, while there is a 
less significant overdensity ~ 380 kpc to the SW. There 
is also an overdensity ~ 1.7Mpc just west of south in the 
outskirts of RXJ1720. 

3.2.3. Substructure in the local kinematics 

Having identified evidence for substructure in both the 
velocity and spatial distributions, we now search for cor- 
relations between the local kinematics and the projected 
local overdensities in the galaxy distribution. As a first 
attempt at this, we utilize the K-test which identifies 
kinematic substructures by comparing the local veloc- 
ity distribution of the N = ^/ji mcm nearest neighbors 
around each galaxy to the global velocity distribution 
(jColless fc Dunnlll996h . The comparison is quantified 
using the Kolmogorov-Smirnov (KS) test, which returns 
the likelihood, Pks, that the local and global velocity dis- 
tributions are drawn from the same parent distribution. 
The overall level of substructure within the cluster can 
be quantified by summing the individual — log[P^s] to 
obtain the re-statistic. The significance of the re-statistic 
is obtained by comparing the observed value to the distri- 
bution of 10,000 remeasurements of the re-statistic which 
have been obtained from samples where the member ve- 
locities are randomly redistributed amongst the posi- 
tions, which are held fixed, thus erasing coherent dynam- 
ical substructure. For A2142 and RXJ1720 we measure 
K-statistics of 836 and 256 at a significance of 5.6cr and 
2.5<7, respectively, indicating that A2142 harbors signif- 
icant dynamical substructure, while RXJ1720 shows ev- 
idence for dynamical substructure, although at a lower 
significance level. 

The results of the re-test for the full field of view for 
A2142 and RXJ1720 are shown in the top right pan- 
els of Figures and HI respectitvely. At the position 
of each member galaxy, we plot a circle with radius 
R oc —log[PKs]- Each circle is color coded so that mem- 
bers with positive i>p CC are red while members with neg- 
ative v pec are blue. Furthermore, those members having 
—log[PKs] values occurring in only 5% of the random 
realizations are plotted with thick lines. Clusterings of 
these large, emboldened circles indicate that there exists 
localized velocity substructure. For comparison to the 
substructures revealed in Section f3. 2. 2 [ we also overplot 
the corresponding surface density contours. 

Within a 1 Mpc radius of the center of A2142 (left panel 
of Figure [§]), there is one region of large clustered circles 
~ 1 Mpc to the NW which is coincident with an exten- 
sion in the galaxy surface density contours. There is 
no significant clustering of large circles coincident with 



the local galaxy overdensity associated with the second 
BCG ~ 180 kpc to the northwest. We note, however, 
that this does not rule out the existence of a compact 
dynamical substructure here, particularly given the en- 
hanced projection effects due to the close proximity to 
the main cluster core which serve to dilute differences in 
the velocity distribution of a real substructure and the 
main cluster. In the cluster outskirts (top right panel of 
Figure [7]), there are three concentrations of large circles, 
two of which coincide approximately with local overden- 
sities in the spatial distribution of galaxies. The most 
significant of these is located ~ 2.6 Mpc to the northwest 
and contains a large number of members with negative 
Up CC . The second lies ~ 2.4 Mpc to the SE, while the 
third, which is not associated with a local overdensity in 
the galaxy distribution, lies ~ 2.3 Mpc due west. There 
is a fourth region ~ 2.4 Mpc to the SW and coinciding 
with a local overdensity in the galaxy distribution which 
contains a more mild concentration of large circles. Fur- 
ther investigation of this region is required to confirm the 
existence of dynamical substructure here. 

Considering the central 1 Mpc region of RXJ1720 
(right panel in Figure [9]) it can be seen that there is 
a clustering of large circles associated with the promi- 
nent secondary galaxy surface density substructure lo- 
cated ~ 550 kpc to the north of the cluster center. There 
is also a smaller clustering of large circles ~ 380 kpc to 
the SW which is coincident with one of the more mi- 
nor substructures in the galaxy surface density distribu- 
tion, indicating there is dynamical substructure there. 
At larger radii (top right panel of Figure H]), we note 
that there are a number of larger circles associated with 
the galaxy overdensity ~ 1.7 Mpc just west of south, al- 
though only a fraction of those are deemed to be sig- 
nificantly larger than expected, indicating that further 
investigation is required before concluding that there is 
local dynamical substructure here. 

The re-test is an ideal tool for detection of substructure 
but provides little information on kinematic properties 
which are key to understanding and characterizing sub- 
structure. For the purpose of characterizing these sub- 
structures, we produce maps tracing the mean velocity 
and velocity dispersion fields which are presented in the 
lower panels of Figures and M for A2142 and RXJ1720, 
respectively. Briefly, these maps were produced by gen- 
erating a grid of 50 x 50 kpc pixels and determining at 
each pixel the weighted mean, u^, 

N 

^= l -^-Jj (!) 

i=l 

and dispersion, a(v w ), 

N 

°{v w ? = — jf (2) 

i=l 

of the N = y/n mcm nearest neighbors, where u pec ,i is 
the peculiar velocity of the ith near neighbor, Wi = 1 — 
(ri/R), Ti is the radial distance of the ith nearest member 
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Figure 7. Top left panel: Isopleths showing the galaxy surface density distribution. The levels are not spaced equally, but are spaced at 
intervals of 5, 10, 20 and 50 galaxies/Mpc 2 between the levels 20-50, 50-100, 100-200 and above 200 galaxies/Mpc 2 , respectively. The 20, 
50, 100 and 200 galaxies/Mpc 2 levels are colored yellow, red, light blue and dark blue, respectively. Green diamonds show the member 
positions. Top right panel: Results of the ft-test where the circle size gives an indication of the difference of local velocity distribution 
compared to the global cluster velocity distribution. Clusterings of large emboldened circles indicate significant departures. Blue and 
red circles are centered on galaxies which have negative and positive peculiar velocities, respectively. Bottom panels: Maps of the mean 
velocity (left) and velocity dispersion (right) fields generated using Equations [l] and [2] The units of the values on the colorbars are kins" 1 . 
Contours in the top right and bottom panels are the same as those shown in the top right panel. 



from the pixel center and R = ri = j\r- 

Examining the results for A2142 in the bottom two 
panels of Figure [7] we note that there is a general co- 
incidence of regions which contain clustered large cir- 
cles seen in the K-test bubble plot and regions where 
either one or both of the and cr(v w ) fluctuate away 
from the overall cluster values. The region of clustered 
large circles ~ 1 Mpc to the NW coincides with a re- 
gion where ~ 750 km s -1 while the most significant 
clustering of large bubbles located ~ 2.6 Mpc to the NW 
coincides with a region where ~ — 600 km s -1 and 



<r(v w ) ~ 400 km s -1 . The substructure ~ 2.4Mpc to 
the SE has mainly negative uHJ values which peak at 
~ — 550 km s -1 just to the north of the local peak 
in the galaxy density. Similarly, the velocity field sur- 
rounding the local galaxy overdensity ^2.5 Mpc to the 
southwest has values of ~ — 450 km s -1 , while the 
substructure ~ 2.5 Mpc due west has ~ 500 kms -1 . 
These maps also reveal that the kinematics immedi- 
ately surrounding the second BCG ~ 180 kpc NW of the 
cluster center are different, with ~ 500 km s -1 and 
a(v w ) ~ 1200 kms" 1 . 
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Figure 8. The same as Figure[7]but for RXJ1720. Here, the lowest isopleth is 10 galaxies/Mpc 



For RXJ1720 (Figure [3]) we find a general agreement 
between the regions which have fluctuations in either one 
or both of and a{v w ) and regions containing the most 
significant local departures from the global velocity dis- 
tribution according to the re-test. The region surrounding 
the local projected galaxy overdensity located ~ 550 kpc 
north of the cluster center shows very complex velocity 
structure. The regions just SW and NE of this substruc- 
ture have ~ — 800 kms -1 and ~ 250 kms" 1 , respec- 
tively, while the the region immediately to the SE of this 
substructure has a(v w ) ~ 1200 kms -1 . The velocity field 
coincident with the substructure ~ 400 kpc to the SW 
has ~ 450 km s , while the velocity field surrounding 
the local galaxy overdensity ~1.7Mpc to the south has 
~ 400 kms -1 , and shows a peak of v w ~ 900 kms -1 
just to the east of the local peak in the galaxy density. 



The re-test and the velocity field maps described above 
are ideal for identifying local kinematic substructure and 
go some way towards classifying the dynamical prop- 
erties of the substructure. However, the data at hand 
allow us to go one step further in the process of vi- 
sualizing and characterizing substructure by producing 
"tomograms" of the galaxy density in peculiar velocity 
slices. These tomograms are shown in Figures [10] and 
[TT1 for A2142 and RXJ1720, respectively. We produce 
maps at nine central velocities where the x — y — v pcc 
distribution has been smoothed with a 3D Gaussian Ker- 
nel with an adaptively varying a xy (the same kernel de- 
scribed in Section f3.2.2[) and a v = 300 kms -1 . The cen- 
tral velocities used in the smoothing are listed in the top 
left corners of Figures [TU] and QT] which show contours 
of the 3D galaxy density in units of galaxies per Mpc 2 
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Figure 9. Zoomed version of the top right panels in Figures [7] and [8] showing the central 1.5 Mpc regions. 



per 1000 km s -1 . The central core is the dominant fea- 
ture in all of the tomograms for RXJ 1720 (Figure ITT]) 
which confirm the complex nature of the dynamics sur- 
rounding the local galaxy overdensity located ~ 550 kpc 
to the N of the core — it appears most significant in the 
v pcc = -1000, -500, 500 and 1000 km s" 1 maps, but is 
relatively insignificant in the v pcc = Okms -1 map. Un- 
like in RXJ 1720, the core does not dominate in all of the 
tomograms for A2142 (Figure IT0"| which are also more 
complicated than their counterparts in RXJ1720. Com- 
parison of the tomograms with those maps presented in 
Figure [7]reveals that those substructures seen in the clus- 
ter outskirts generally appear to be fairly well isolated 
structures in velocity space, too. The tomograms show 
that the conglomeration of significant circles ~ 1 Mpc 
to the NW of the core appears as a feature in the 
^pcc > 500 km s -1 maps, while the local overdensity as- 
sociated with the second BCG appears most prominently 
at tJ poc = 1500 and 2000 km s" 1 . 

3.2.4. Characterizing the substructures with KMM 

Having identified substructure with a number of tech- 
niques in Sections 13.2.21 and 13.2.31 we are now in a posi- 
tion to characterize the dynamical properties of the sub- 
structures. To achi eve this, we use the Ka ye's M ixture 
Model algorithm of lAshman et all (KMM: ]l994h to in- 
vestigate the dynamics within the regions shown in Fig- 
ures [T2] and [T3J The siz es of these regions were chosen 
to contain ~ l-5\/Ar moin = 50 and 30 nearest neighbors 
for A2142 and RXJ1720, respectively. For each of these 
regions, the KMM algorithm was used to fit a minimum 
of two Gaussians to the velocity distribution. The first 
Gaussian component models the main cluster velocity 
distribution while the remaining Gaussian components 
model the superimposed substructure velocity distribu- 
tion. The KMM algorithm requires initial estimates of 
the mean (//), variance (a 2 ) and the fraction of galax- 
ies (/) belonging to each Gaussian velocity component. 
The initial estimate for the main cluster component is set 
to (/ii,in, ci,in, /i,i n ) = (0 km s" 1 , 1000 km s _1 , 0.8) , while 



the initial estimates for the remaining Gaussian compo- 
nents are determined by inspection of the velocity distri- 
bution and are listed in Tables [5] and El We determine the 
significance of each bimodal fit by employing a paramet- 
ric Bootstrap method (see lOwers et al.l I201TI ). Briefly, 
we produce 1000 random resamplings of a Gaussian with 
50 (for A2142) or 30 (for RXJ1720) data points and \i 
and a set to the values derived from the entire cluster 
sample. The KMM algorithm is then used to fit two 
Gaussians to the resampled distribution with the initial 
inputs set to the outputted KMM fits to the observed 
data. A measure of the improvement of the fit in go- 
ing from a unimodal to bimodal Gaussian distribution is 
given by the Likelih ood Ratio Test Statistic (LRTS; see 
lAshman et"aTJ I1994T ). The distribution of LRTSs pro- 
duced by the resampled distributions can be compared 
to the observed LRTS to determine if the results of the 
fit to the observed distribution can be caused by random 
fluctuations in the data and this is listed as a P-value 
in Table [5] and El In a number of cases (e.g., S5 and 
S6 for A2142 and SI for RXJ1720) a three-mode par- 
tition was also fitted to the data. For these cases, the 
P-value is evaluated by comparison of the trimodal fit to 
the bimodal one. 

3.3. Summary of the substructure properties 

Having identified many significant substructures in 
both A2142 and RXJ1720, and characterized them dy- 
namically using the KMM method, we now briefly discuss 
each of them and, in particular, their properties that are 
most pertinent to our study. 

3.3.1. A2142 

• SI: This substructure is coincident with the sec- 
ond BCG ~ 180 kpc to the northwest of the cen- 
tral BCG. There is a local overdensity in the pro- 
jected galaxy distribution here, although the re- 
test does not reveal strong evidence that the lo- 
cal velocity distribution differs significantly from 
the global one. The substructure is prominent as 
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Figure 10. Tomograms for A2142 revealing overdensities in position-velocity space centered at the velocities listed in the upper left of 
each panel. The spatial smoothing uses the same adaptive kernel as that used in the top left panel of Figure[7] while the velocity smoothing 
uses a fixed sigma of 300 km s" 1 . The green diamonds show the positions of galaxies lying within ±600 km s 1 of the central velocity, v pec . 
The contours show the galaxy density and start at 5 galaxies/Mpc 2 /(1000 km s _1 ) with an interval of 5 galaxies/Mpc 2 /(1000 km s — 
up to the 50 galaxies/Mpc 2 /(1000 kms" 1 ) level (red contour). Above 50 galaxies/Mpc 2 /(1000 km s -1 ), for clarity, the interval is 
10 galaxies/Mpc 2 /(1000kms- 1 )- 



Table 5 

KMM fits to the regions shown in Figure IT2l for A2142. 



Region 


(X, Y) 


Substructure Input 


Main output 


Substructure output 


P-value 




(kpc, kpc) 


(M, a, f) 


(/*■ *, f) 


(M, t, f) 





51 (150, 90) (1800, 300, 0.2) (-232.8, 938.7, 0.86) (1732.5, 224.2, 0.14) 0.434 

52 (600, 763) (1500, 300, 0.2) (186.3, 638.3, 0.82) (1681.0, 368.8, 0.18) 0.007 

53 (2007, 1567) (-600, 400, 0.2) (500.1. 688.7, 0.18) (-450.8, 442.5, 0.82) < 0.001 

54 (-1255, -2038) (-900, 300, 0.2) (114.7, 1047.5, 0.54) (-913.4, 178.5, 0.46) < 0.001 

55 (2072, -1494) (-1500, 300, 0.2) (272.0, 649.9, 0.72) (-1608.3, 153.1, 0.28) < 0.001 

55 (2072, -1494) (-1500. 300, 0.2), (600, 200, 0.1) (24.2, 765.6, 0.44) (-1620.4, 144.1, 0.26), (658.8, 68.9, 0.2) 0.007 

56 (2327, -180) (1200, 300, 0.2) (-701.1, 522.5, 0.40) (1156.5, 357.8, 0.60) < 0.001 

56 (2327, -180) (1200, 300, 0.2), ( -400, 200, 0.1) (-1212.4, 423.8, 0.14) (1147.8, 362.8, 0.60), (-433.7, 220.7, 0.26) 0.32 

57 (-1725, -100) (900, 300, 0.2) (-185.0, 815.2. 0.78) (902.6, 172.2, 0.22) 0.11 
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Table 6 

KMM fits to the regions shown in Figure [13] for RXJ1720. 



Region 


(X, Y) 


Substructure Input 


Main output 


Substructure output 


P-value 




(kpc, kpc) 


a , f) 


(M, t, f) 


Cf*, «t, /) 




SI 


(-50, 576) 


(700, 300, 0.2) 


(-849.7, 594.4, 0.67) 


(736.0, 184.2, 0.33) 


0.008 


SI 


(-50, 576) 


(700, 300, 0.2), (-700, 300, 0.2) 


(-906.0, 659.9, 0.50) (733.8 


185.6, 0.33), (-658.2, 127.9, 0.17) 


0.72 


S2 


(443, -144) 


(-900, 300, 0.2) 


(-42.4, 1145.8, 0.77) 


(-851.0, 99.2, 0.23) 


0.082 


S3 


(302, -1706) 


(800, 300, 0.2) 


(-819.8, 404.8, 0.23) 


(818.4, 512.5, 0.76) 


0.009 


S4 


(-1225, 181) 


(-400, 300, 0.2) 




(-291.1, 577.3, 1.00) 
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Fig ure 12. Top panel: The black circles show regions selected for the KMM analyses based on substructure revealed in Figure [71 and 
1101 Black contours show galaxy surface density as in Figure [jj The stars are color coded to match the corresponding velocity components 
shown in the lower panels. Black points reveal the positions of the remaining cluster members. 
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Figure 13. As for Figure \U\ but for RXJ1720. 
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an overdensity in the T7 pcc = 1500 and 2000 kms -1 
tomograms, consistent with a grouping of galax- 
ies around the velocity of the second BCG (u P oc = 
1760 kms -1 ). The KMM analysis returns a mean 
Wpoc = 1732.5 kms -1 and dispersion <7 vpec = 
224.2 kms -1 for a group of ~ 7 — 10 galaxies sur- 
rounding this substructure, although a bimodal fit 
is not favored. The lack of evidence for dynamical 
substructure in the re- and KMM tests is proba- 
bly due to Si's proximity to the cluster core where 
projection effects caused by the main cluster pop- 
ulation are large. This is likely the remnant core of 
a more massive system which has had its less well 
bound members stripped during its infall. 

S2: This is a more loose agglomeration of spi- 
rals. There doesn't appear to be a dominant early 
type member within this substructure and there 
is only a mild excess in the galaxy surface den- 
sity here. The velocity field is significantly posi- 
tive in this region and the tomograms reveal lo- 
cal overdensities here which are most prominent at 
TJpcc = 1000 and 1500 km s -1 . The KMM fit reveals 
a significant bimodality where the substructure has 
mean velocity v pcc — 1681.0 kms -1 and dispersion 
of cr vpcc = 368.8 km s" 1 . lOkabe fc Umetsul pOOl 
find a mildly significant peak in their projected 
mass maps at approximately the same position. 

S3: This substructure is clearly distinguished by 
the re-test. It lies ~ 2.5 Mpc from the cluster cen- 
ter and it is conspicuous as a local overdensity in 
the projected galaxy distribution, as features in the 
velocity field maps and as an overdensity in the 
Vp C c = — 500 kms -1 tomogram. The KMM analy- 
sis returns significant bimodality and the substruc- 
ture has mean velocity v pcc = —450.8 kms -1 and 
dispersion (T vpoc = 442.5 kms -1 , consistent with a 
large group or small cluster. 

S4: Another substructure in the outskirts located 
~ 2.4 Mpc from the center. This substructure man- 
ifests itself as a local overdensity in the projected 
galaxy distribution and in the zJ pcc = — 1000 kms 
tomogram. The KMM test reveals significant bi- 
modality and the substructure has mean veloc- 
ity f poc = —913.4 kms -1 with dispersion a = 
178.5 km s , consistent with a group-sized sub- 
structure. 

S5: This substructure is a significant local overden- 
sity in the projected galaxy distribution although 
the re-test does not show conclusive results for local 
velocity substructure. The substructure is seen as 
an overdensity in the v pcc = — 1500 and 500 kms -1 
tomograms. The KMM analysis reveals the ve- 
locity distribution is trimodal with two low dis- 
persion substructures with mean velocity v pcc — 
— 1620.4and658.8kms -1 and dispersions cr v = 
144.1 and 68.9 kms -1 , respectively. The low veloc- 
ity dispersions indicate that the substructures have 
low mass and their large peculiar velocity separa- 
tion indicates they are not associated. The projec- 
tion effects caused by these two separate structures 



aligned along the line of sight have enhanced the 
projected galaxy density here. 

• S6: There is only a mild increase in the local 
galaxy surface density in this region, however the 
re-test shows that the local kinematics are signif- 
icantly different from the global kinematics. The 
KMM analysis showed that the reason for this 
difference was a strong bimodality with two sub- 
structures separated by a peculiar velocity of ~ 
1860 kms -1 . The two substructures have mean ve- 
locity Vpcc = —701, 1 and 1156.5 kms -1 and disper- 
sion CT V pcc = 522.5 and 357.8 kms -1 , respectively. 
The tomograms show that the negative peculiar 
velocity component is fairly dispersed, while the 
Vpcc ~ 1000 kms -1 component does appear to be 
more concentrated in this region, although it is not 
as prominent as other substructures such as S3. 
The galaxies in this region may be part of coher- 
ent substructures which are infalling from the sur- 
rounding large scale structure. 

• S7: This substructure is prominent as a local over- 
density in the projected galaxy distribution, al- 
though the re-test does not reveal any significant 
velocity substructure here. There are local den- 
sity peaks here in the v pcc = 500 and 1000 kms -1 
tomograms (Figure [T0|) and the KMM algorithm 
fits a substructure with mean velocity v pcc = 
902.6 kms -1 and dispersion cr vp cc = 172.2, al- 
though bimodality is not strongly favored. The low 
dispersion indicates a low mass, although given the 
minimal evidence for velocity substructure, the lo- 
cal peak in the galaxy density here may simply due 
to projection effects from unrelated galaxies. 

3.3.2. RXJ1720 

• SI: This substructure is a well defined lo- 
cal peak in the galaxy surface density distribu- 
tion although it has unusual kinematic proper- 
ties. The re-test reveals significant local kine- 
matic substructure in this region, while the to- 
mograms reveal local overdensities in this region 
at vpcc = -500, -1000, 500 and 1000 kms -1 but 
not at tJpcc = Okms -1 . The KMM analysis re- 
vealed a bimodal fit was preferred over uni- and tri- 
modal fits with the two components having mean 
velocity v pcc = —849.7 and 736.0 kms -1 and dis- 
persions of cr V pcc = 594.4 and 184.2 kms -1 . The 
second brightest member in the core of the clus- 
ter ha s i> Dec = 933 kms -1 and is within the SI re- 
gion. lOkabe et al.l (|2010l) find a second peak in 
their weak lensing convergence maps which is co- 
incident with this substructure, providing further 
evidence that this substructure is real. We tenta- 
tively point to a mild excess in the X-ray residuals 
in the right panel of Figure Q21 which is coincident 
with this structure. 

• S2: This substructure has a mild enhancement in 
the projected galaxy distribution and a local ve- 
locity distribution which differs from the global one 
according to the re-test. The KMM analysis reveals 
that this difference is caused by an extremely low 
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velocity dispersion (<7 vpcc = 99.2 km s 1 ) substruc- 
ture at a velocity u poc = —851.0 km s _1 , although 
the improvement in the fit in going from a uni- 
modal to bimodal distribution is marginal. Given 
the small velocity dispersion and mild enhance- 
ment in galaxy surface density, this substructure 
is probab ly dynamical l y insi gnificant although we 
note that lOkabe et al.l (|2010f ) find a weak enhance- 
ment at roughly this position in their lensing maps. 

• S3: This substructure shows a mild local peak in 
the projected galaxy distribution and the velocity 
field shows a preference for positive peculiar veloc- 
ities here. While the tomograms do not reveal sig- 
nificant local peaks for the maps centered around 
positive peculiar velocities, the KMM analysis re- 
veals that the velocity distribution here is bimodal 
with the majority of the galaxies allocated to a 
component centered at u poc = 818.4 km s _1 with 
dispersion <7 vpoc = 512.5. This large dispersion, 
if real, indicates a significant mass in this region. 
However, visual inspection of the velocity distribu- 
tion in the lower left panel of Figure Q2] reveals that 
there may be some contamination from galaxies at 
v P cc ~ 1500 km s -1 which may be inflating the dis- 
persion. More complete spectroscopy is required in 
this region to better disentangle the velocity com- 
ponents and confirm the large dispersion measured 
here. 

• S4: The K-test revealed a number of large bub- 
bles in this region which indicates evidence for 
velocity substructure. There is no strong galaxy 
surface density enhancement, although the tomo- 
gram centered at u pcc = — 500kms _1 shows a lo- 
cal peak here while the velocity dispersion in this 
region appears to be low. However, the KMM 
analysis failed to fit a second velocity substructure 
to the data in this region, although the best fit- 
ting single Gaussian does have a low velocity dis- 
persion (<7 V pec = 577.3 km s _1 ) and is centered at 
v P cc = — 291. lkm s -1 indicating a departure from 
the global velocity distribution here. More data 
are required in this region in order to confirm the 
existence of real dynamical substructure. 

4. DISCUSSION 

The principal objective of this study is to find obser- 
vational evidence for merger related substructure in the 
sloshing cold front clusters A2142 and RXJ1720. Fur- 
thermore, we wish to determine if the cold fronts ob- 
served in the cores of these clusters can be related to 
merger activity, therefore confirming that cold fronts 
are excellent signatures of merger activity. The re- 
sults of our substructure tests reveal that A2142 and 
RXJ1720 harbor a number of group-scale substructures 
which manifest themselves as enhancements in the lo- 
cal galaxy surface density (Section l3.2.2p . and as regions 
of localized velocity substructure (Section 13 . 2 .3[) . The 
existence of these group-scale substructures is consis- 
tent with a minor merger origin of the cold fr o nts, a s 
seen in the simulations of iTittlev fc Hcnrikscn (2005), 
Ascasibar fc Markevitchl (|2006t see also iRoediger et al.l 
20111 for a more specific case) where a low mass group 



perturbs the main cluster core and induces cold fronts. 
The problem now, particularly for A2142, is to deter- 
mine which of the many substructures detected is the 
most likely perturber. As a first step, we note that sub- 
structures with large i? pro j are unlikely to have made a 
pericentric passage and, thus, are unlikely to be the per- 
turber. We outline the evidence for this statement in 
subsection 14.11 and discuss the substructures which we 
consider to be the most likely perturbers in A2142 (SI 
and S2; subse ction EZEJ) and RXJ1720 (SI and S3; sub- 
section E221) ■ 

4.1. Ruling out the substructures least likely to be 
perturbers 

First, we consider the apocentric distance of a sub- 
cluster after its initial pericentric passage. We obtain 
an esti mate for this di s tance by considering the simula- 
tions of iTormen et al.l (|2004f) where the average apoc- 
entric distances for an ensemble of simulated clusters 
with merging subhalos of varying masses is shown in 
their Figure 11. We only c onsider results for th e "self 
bound" satellite particles in ITormen et al.l (|2004| ). since 
these results are more likely to reflect our observed sub- 
structures. Two important features can be noted here: 
(i) On average, the subhalos do not travel farther than 
the main cluster's virial radius, Rv, after pericenter 
and (ii) Due to dynamical friction, larger mass sub- 
halos have apocentric distances as little as ~ 0.5i?y- 
Thus, it is highly unlikely that any substructure ob- 
served at radii larger than the cluster virial radius has 
passed pericenter — they are far more likely to be groups 
which are falling onto the cluster for the first time. For 
A2142 and RXJ1720, we assume that the virial radius 
is i?vir — -R200 = 0-17<t v /H(z) where i?2oo is the radius 
where the density of the cluster is 200 times the critical 
density at redshift z, o~ v is the cluster velocity disper- 
sion and H(z) = H ^/£l m {l + z) 3 + fly (]Carlberg et al.l 
119971) . This results in i? vir = 2.36 and 2.02 Mpc for A2142 
and RXJ1720, respectively, and we set this as an initial 
constraint on the i? pro j within which the perturber must 
lie. These limits are plotted as cyan circles in Figures IT4l 
and [151 For RXJ1720, these limits do not strongly rule 
out any of the substructures as potential perturbers, al- 
though S3 lies close to the virial radius in projection. For 
A2142, the substructures S3 and S5 are outside of the 
Rv and are therefore ruled out as potential perturbers, 
while S4 and S6 lie at approximately the virial radius. 
Those substructures at the virial radius in A2142 (S4, 
S6) are unlikely to be the perturbers because (i) Projec- 
tion effects mean that the observed radius gives the min- 
imum distance from the cluster center — if the substruc- 
ture lies at some distance in front or behind the cluster, 
then it is likely to lie outside the virial radius and (ii) 
The limit imposed here is an upper limit on the radius. 
More massive systems which have passed close to the 
cluster core are stro ngly affected by dynamical friction 
(|Tormen et al.l [20041 ) and thus have smaller apocentric 
distances. 

Further evidence that these substructures at large 
-ftproj (<; Rv) are unlikely to have passed pericenter 
comes from considering the timescales required to pro- 
duce the cold front pairs seen in the Chandra images 
(Figures [14] and [15]) . Given that the fluid velocities in- 
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volved are small compared t o the Kepler speed (e.g., 
IAscasibar~ Sz Markevitch 200(f), the sloshing motions that 
give rise to cold fronts are likely to be reasonably well 
described as superpositions o f internal gravity wav es ex- 
cited by merging subclusters (jChurazov eTalBonal) . Lo- 
cally, the oscillation frequency of these waves is the 
Brunt- Vaisala frequency 

/ g 3 d In s 13 d\ns 

V r 5 d In r V 5 din r ' 

where g = GM (r) / r 2 is the acceleration due to grav- 
ity, r is the radius, s is the entropy index and wk is the 
Keplerian orbital frequency. Thus ubv is expected to 
provide an estimate of sloshing frequencies. The radial 
entropy profi les for both A2142 and RX J 1720 have been 
measured by iCavagnolo et al.l (|2009l ) and they have log- 
arithmic slopes of ~ 1.22 and 1.39 at the positions of the 
cold fronts, respectively Therefore, the approximation 
ojbv — — Vcitc/r is s ufficient for our purpose s. This 
justifies the estimate of ISimionescu "etall (pOlOl ). which 
employs the radial dependence of the frequency to de- 
termine the time required for an adjacent pair of cold 
fronts to get out of phase by tt radians, i.e., to appear 
on opposite sides of the cluster center. Even in the lin- 
ear approximation, a more thorough treatment requires 
determination of the spectrum and phasing of the inter- 
nal modes excited by a merging subhalo, then following 
the time development of this superposition of internal 
gravity waves. Apart from nonlinear effects, this leaves 
considerable uncertainty in our time estimates. 

To determine uif., we measure the cold front radii from 
the central brightest cluster galaxy which is assumed 
to mark the potential minimum and assume a singular 
isothermal sphere (SIS) density distribution meaning the 
circular velocity is simply V c i rc — V2a v . We also as- 
sume that the time required for the cold front pairs to 
become tt radians out of phase in their orbits is equiv- 
alent to the time since the perturber's pericentric pas- 
sage. For A2142, the inner and outer cold fronts have 
radii -R pro j,m = 87kpc and -R pro j,out = 326 kpc, respec- 
tively, and V c i IC = 1407 km s -1 . The time required for 
the cold fronts to become out of phase by tt radians is 
t = 7T/(u>k,in — Wfe )OU t) — 0.26 Gyr. The same calcula- 
tion for RXJ1720 with J? pro j,m = 69 kpc, i? pro j,out = 
185 kpc and Km, = 1247 km s -1 gives r ~ 0.27 Gyr. 
IRoediger et al.l (120111) found that t he cold front age cal- 
culated by ISimionescu "etaLl (l2010l) was underestimated 
by a factor of 3.5 when compared to their simulations. 
We have also c alculated r for several of the co ld fronts 
in Figure 7 of lAscasibar fc Markevitch! (|2006ft . finding 
that this method generally underestimates the time since 
pericentric passage by fact ors of 3-4, consistent with the 
findings of IRoediger et all . Therefore, we estimate that 
the "true" time since pericentric passage for both A2142 
and RXJ1720 is r truc ~ 0.8 - l.lGyrs. We note that 
viewing angles and uncertainties in the phasing of the 
waves mean that these estimates give only a rough guide 
for the time since pericenter. Furthermore, we note that 
the cold fronts that we observe in A2142 and RXJ1720 
occur on much larger scales than the cold fronts seen at 
similar times in the simulations. This may reflect differ- 
ences in the physical properties of A2142 and RXJ1720 
compared to the simulated clusters, e.g., the shape of 



the underlying potential will affect the oscillation fre- 
quencies. The spatial and temporal scales of the fronts 
are also likely to be affected by properties of the per- 
turber, including its mass and impact parameter. Nev- 
ertheless, the timescales derived above provide useful es- 
timates for our purposes and prov ide strong support fo r 
the radial constraints derived from lTormen et all (|2004|) . 
For a substructure to have reached the virial radius af- 
ter pericenter on these timescales, it would need to have 
traveled at an average radial velocity since pericenter of 
- 2100 - 2890 and 1790 - 2460 km s" 1 for A2142 and 
RXJ1720. While velocities of this order are expected 
at pericenter during a merger, they do not persist on 
gigayear timescales since the subcluster velocity quickly 
decays after pericenter due to the effects of gravity and 
dynamical friction as it travels outward. This places 
strong constraints on whether substructures S4 and S6 
in A2142 have had time to have passed pericenter and 
traveled to their current i? pro j. Therefore, we disregard 
them as being potential perturbers. 

Having ruled out S3, S4, S5 and S6 in A2142, we now 
consider substructure S7 which we rule out for the fol- 
lowing reasons. First, while it lies within the virial ra- 
dius, it has a relatively large radius (-ff pro j ~ 1.7 Mpc) 
which, as noted above for S4 and S6, gives its minimum 
distance from the cluster center. Furthermore, as noted 
in Section 13.31 while there is a local excess in the pro- 
jected galaxy density, there is little evidence for any cor- 
responding localized velocity substructure. This means 
that if there is a real substructure in this region, it is 
probably not very massive. A further possibility is that 
the excess in galaxy surface density is due to projection 
of galaxies which are not physically associated and there- 
fore that S7 is not a bona fide substructure. Given its 
uncertain nature and large i? pro j > we find it unlikely that 
S7 has passed pericenter and, if it has, it is probably 
not massive enough to have caused the perturbation re- 
quired to produce the cold fronts with such large radii 
and contrast. Therefore, we also dismiss S7 as a potential 
perturber. 

The radial limits derived above for RXJ1720 placed no 
strong limits on any of the substructures shown in Fig- 
ure Uni However, we rule out the substructures S2 and 
S4 because they are unlikely to be significant substruc- 
tures, if real associations at all. Starting with S2, this 
substructure consists of a grouping of a small number 
of galaxies with similar velocities and a small dispersion 
of only ~ 100 kins" 1 . There appears to be a substruc- 
ture at approximat ely this location in the lensing maps of 
IQkabe et al.l ((2010), although it is not concentrated and 
is only detected at the 3cr level. The presence of sub- 
structure detected using our spectroscopica lly confirmed 
memb ers and also in the lensing maps of IQkabe et al.1 
(2010) provide strong support that S2 is a real substruc- 
ture. However, the low velocity dispersion and low sig- 
nificance in the lensing maps imply that S2 has low mass 
and is therefore not likely to have played an important 
role in perturbing the cluster core. Considering the S4 
substructure, while there are hints of localized velocity 
substructure in this region according to the K-test and 
the velocity fields, there is only a weak excess in galaxy 
surface density here. Furthermore, the KMM analysis 
failed to identify any bimodality in this region. Given 
the uncertainty in the existence of velocity substructure 
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and the lack of a compact 2D substructure, it is unlikely 
that S4 is a real substructure. Deeper, more complete 
spectroscopic coverage of RXJ1720 would be required to 
determine if there is a real substructure there. 

4.2. Most likely perturbers 
4.2.1. A2142 

We begin our discussion with the SI substructure sur- 
rounding the second BCG at ~ 180 kpc to the north- 
west of the first ranked BCG. The KMM fit returns 
v P cc — 1732.5 and a v — 224.2kms _1 , although a bi- 
modal fit to the local velocity distribution was not fa- 
vored statistically over a unimodal one. However, the 
surrounding local overdensity in the galaxy surface den- 
sity, along with the compact local overdensity in the 1500 
and 2000kms _1 tomograms (Figure ITO)) indicate that 
there is a kinematically compact grouping of ~ 7—10 
galaxies associated with this BCG. Any kinematic dis- 
tinction is likely diluted by the projection of the main 
cluster core galaxies onto SI due to its close proximity in 
projection to the core. BCGs are expected to be found at 
the centers of groups and clusters thus the existence of a 
second BCG provides strong evidence that SI is a merger 
remnant. Given the relatively high (although not ex- 
treme) positive peculiar velocity of SI, we speculate that 
the majority of its motion is aligned with our line of sight. 
In the scenario where SI is the perturber responsible for 
the cold fronts, it must have already passed pericenter. 
If the estimated time since pericenter of ~ 1 Gyr derived 
above is correct, then SI is unlikely to have passed apoc- 
enter and therefore we surmise that it is currently located 
on the far side of the cluster on its way to apocentcr. 

Further evidence for a perturber which has an orbit 
which is highly inclined to the plane of the sky, such 
as that proposed for SI, comes from the morphology 
of the cold fronts. The Chandra data for A2142 reveal 
two prominent cold fronts (|Markevitch et al. 2000) strad- 
dling the cool core (Figure [14]) . There is no clear spiral 
structure seen in the residual map (right panel of Fig- 
ure I14[) , n or in the temperature maps available in th e 
literature (|Markevitch et al.1 120001 lOwers et al.l l2009bf l. 
although there are fairly linear structures joining the in- 
nermost cold front to the outermost one. The morphol- 
ogy of the cold fronts is consistent with a spiral struc- 
ture viewed with its axis almost perpendicular to our 
line of sig ht, e.g., as shown in the middle panel of Fig- 
ure 19 in lAscasibar fc Markevitchl (|2006T l . In the sim- 
ulations, the spiral structure is formed when angular 
momentum is transfered to the displaced core gas from 
a perturber which has underg one an offcenter passage 
(Ascasibar & Markcvitch] 120061 ). Because the angular 
momentum responsible for the spiral structure originates 
with the offcenter passage of the perturber, the spiral and 
perturber orbit must be coplanar. Based on this and the 
evidence that we are viewing a spiral structure with its 
axis almost perpendicular to our line of sight, we surmise 
that the perturber is orbiting in a plane which is highly 
inclined to the plane of the sky. 

We also note here that the central BCG has a peculiar 
velocity of v pcc = 212 ± 65kms~ 1 with respect to the 
cluster velocity distribution. This difference is signifi- 
ca nt at the 2.9a l evel a ccording to the method outlined 
in lOegerle fc Hilll (|2001l ). Assuming that the BCG lies at 



the bottom of the cluster potential as traced by the dark 
matter, we interpret this as evi dence of dark matter core 
oscilla tions. As pointed out by [Ascasibar & Markcvitch 
(2006) the observation of BCG peculiar velocities is nat- 
urally encompassed in their gasless minor merger sce- 
nario. In that scenario, a purely gravitational pertur- 
bation excites oscillations in both the dark matter and 
gas cores leading to cold fronts when the gas decouples 
from the dark matter. However, in the case where a 
weak shock or other purely hydrodynam i cal disturbance 
(iChurazov et al.ll2003l iFujita et al.ll2004 iRoediger et al.l 
1 2 1 lh perturbs the core and generates the cold fronts, 
only the gas is expected to oscillate, thus no significant 
BCG peculiar velocity should be observed. We inter- 
pret the BCG peculiar velocity as evidence that the cold 
fronts are generated by the gravitational disturbance of 
a minor merger. Furthermore, the detection of a line of 
sight velocity offset lends support to our assertion that 
the perturber's orbit is inclined to the line of sight. 

As an alternative, we consider the S2 substructure, 
which is a fairly loose group of galaxies at a projected 
distance of ~ 970 kpc from the central BCG, has a 
mean peculiar velocity of v pcc = 1680 km s -1 and dis- 
persion a ~ 370 km s . Supporting evidence for S2 be- 
in g a real substructure com es from the lensing analysis 
of IQkabe fc Umetsul (]2008l ) who detect a signal in their 
projected mass map which is significant at the 3cr level 
and is approximately at the position of S2. Similar to 
SI, the positive peculiar velocity of S2 implies that it 
is currently travelling away from us and, assuming that 
it has not passed apocenter, that its initial trajectory 
began in the foreground to the southeast of the main 
cluster core. It passed pericenter on the SE before being 
deflected behind the cluster core and to the northwest 
where it traveled to its current position. 

4.2.2. RXJ1720 

The most viable perturber candidate for RXJ1720 is 
the SI substructure lying ~ 550 kpc north of the clus- 
ter core. This substructure is conspicuous as an ex- 
cess in the galaxy surface density maps and shows ev- 
idence for localized kinematic substructure. A more de- 
tailed inspection of the velocity distribution surround- 
ing SI in Section 13.2.41 revealed significant bimodality 
with two substructures separated by ~ 1586 km s -1 , nei- 
ther of which is centered at the systemic cluster velocity 
(fix = -849.7kms- 1 and ^2 = 736.0 kms" 1 ). While 
the velocity structure is complicated here and is hard to 
interpret, we are confident that there exis ts a massive 
substr ucture here, particularly given that IQkabe et al.l 
(2010) find significant substructure in their lensing maps 
which is spatially coincident with SI. 

Here, we consider a scenario where SI has perturbed 
the core during its pericenter passage and is currently 
moving outwards towards the north. Considering the 
Chandra residual map generated from the Chandra data 
in the right panel of Figure [T5] a clear spiral struc- 
ture ca n be seen which is also seen in the temperature 
maps o f iMazzotta fc Giacintuccil (|2008l ) and lOwers et al.l 
(2009b|). Following on from the discussion on the spiral 
morphology in Section 14.2.11 this means we must be ob- 
serving the spiral structure which has an axis that is 
aligned roughly with our line of sight. We also note that, 
in contrast to A2142, we see no significant peculiar ve- 
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locity offset in the BCG, consistent with the sloshing 
occurring mainly in a plane perpendicular to our line 
of sight. By extension, the axis of the perturber's or- 
bit is roughly coincident with our line of sight. Fur- 
thermore, following the spiral from the outermost cold 
front to the core, the windin g direction is clockwise. 
iMarkevitch fe Vikhlininl (|2007l) noted that the orbit of 
the perturber and that of the gas generating the spiral 
pattern shou ld exhibit prograde rotation, a s seen in the 
simulations (lAscasibar fc Mark evitchl I2006L iPoole et al.l 
l2006LlRoediger et al.ll2011l ). This implies that the SI sub- 
structure follows a clockwise orbit from our perspective. 
We propose that SI has traveled from the south, passed 
the core region on the eastern side providing the neces- 
sary perturbation for the generation of the cold fronts 
and is now headed to the north. One issue with this in- 
terpretation is the proximity of SI to the cluster core. If 
the time since pericenter of 1 Gyr derived above is cor- 
rect and the perturber's orbital plane is very close to 
the plane of the sky, then SI should have had time to 
travel much further than 550 kpc. Of course, this asser- 
tion relies critically on the velocity of SI which is not well 
constrained. However, even assuming a relatively small 
radial velocity of 1000 km s -1 since core passage implies 
that SI should be ~ IMpc from the cluster core. This 
issue is somewhat relieved if we take the lower value of 
''"true = 0.8 derived above and assume that Si's orbit has 
some component parallel with our line of sight. 

Alternatively, if we have underestimated the time since 
pericenter, then it is possible that SI has passed apocen- 
ter and is currently headed towards its second pericenter. 
This also relieves the tension between the observed scales 
of the cold fronts in RXJ1720, which occur at much larger 
radii when compared with those formed on ~ 1 Gyr 
post-pericenter timescal e s see n in the simulations of 
lAscasibar fc Markevitch l (120061). For example, in the 3. 8 
Gvr panel of Figure 7 in lAscasibar fc Markevitchl (|2006f) . 
i.e., ~ 2.4 Gyr after pericenter and just prior to the sec- 
ond pericenter, the most prominent cold fronts have radii 
~ 80 kpc and ~ 160 kpc, simil ar to the cold fronts ob- 
served in RXJ1720. However, lAscasibar fc Markevitchl 
(|2006f ) note that the contrast in density and temperature 
at the cold fronts decreases with time by a factor of 1.5 at 
3.8 Gyrs, co mpared with the fact or of ~ 2 jump observed 
in RXJ1720 (|Owers et al.ll2009bt) . This difference in con- 
trast may reflect a difference in the mass, the gas content 
and/or the impact parameter of the perturber between 
RXJ1720 and the simul ations presented in Figure 7 of 
Ascasibar & Markevitch. 

The S3 substructure located at ~ 1.7 Mpc south of 
the cluster core may also be considered a potential per- 
turber. The KMM fits revealed that the velocity of 
this substructure is ~ 820kms _1 and it has dispersion 
a r~j 512 kms . This dispersion is fairly large and com- 
parable to that expected of a reasonably large group or 
poor cluster. However, visual inspection of the velocity 
distribution (lower left panel in Figure [T3|) reveals that 
the 820 km s _1 partition has two peaks and thus that this 
partition may be contaminated by galaxies which are not 
associated with the substructure. This contamination 
may be responsible for inflating the velocity dispersion. 
More data are needed to disentangle any contamination 
and to confirm this velocity dispersion. In any case, there 



is clearly a dynamical substructure in this region. In the 
scenario where S3 is the perturber, we propose that it 
has traveled from the north and is currently heading to 
the south. Based on the arguments above on the winding 
direction of the spiral structure seen in the right panel of 
Figure [15j we suggest that S3 passed the cluster core on 
the western side and is currently heading towards apoc- 
enter for the first time. 

5. SUMMARY AND CONCLUSIONS 

We have presented comprehensive MMT/Hectospec 
and Kcck/DEIMOS multi-object spectroscopy of the cold 
front clusters A2142 and RXJ1720. These clusters were 
selected from the "relaxed appearing" s ubsample of cold 
front clusters from lOwers et al.l (|2009bl ) with the specific 
aim of testing the hypothesis that minor mergers cause 
the observed cold fronts. To do this, we used the spa- 
tial and kinematic information contained in samples of 
956 and 400 spectroscopically confirmed cluster members 
for A2142 and RXJ1720, respectively, to search for sub- 
structure. In both of the clusters we find group-sized 
substructure and we identify likely perturbers based on 
comparisons between the observed cold front morpholo- 
gies and those seen in the simulations. Our results are 
consistent with the cold fronts being caused by sloshing 
induced by the gravitational perturbation caused by a 
minor merger. We stress, however, that more detailed 
merger histories and better orbital constraints will likely 
require cluster-specific simulations aimed at reproducing 
both the X-ray morphology and the positions and kine- 
matic properties of the detected substructures. Nonethe- 
less, we conclude that our original stated go al of showing 
cold f ronts are related to merger activity (jOwers et al.1 
I2009d |al|b1) has been met. 
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Figure 14. Left panel: SDSS gri image of A2142 with X-ray contours from ROSAT (cya n) and Chandra (magenta) overlaid. The white 
dashed circles show the same regions used to define the KMM substructures in Figure IT21 The cyan dashed circle shows the cluster virial 
radius (Ry = 2.36 Mpc) beyond which we do not expect to find substructure which has passed pericenter. Right panel: Residuals from an 
exposure and background corrected 0.5-7 keV Chandra image after subtraction of the azimuthal average at each radius. There is no clear 
spiral structure joining the outer cold front to the cluster core in this residual map. 
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Figure 15. Left and right panels: The same as Figure Il4l but for RXJ1720. ROSAT contours are not shown in the left panel, and the 
radius of the cyan dashed circles is Ry = 2.02 Mpc. 
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